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LONG-TERM  GOALS 

The  long-term  scientific  goal  of  my  research  is  to  better  understand  the  distribution  of  phytoplankton 
in  the  world's  oceans  through  remote  sensing  their  influence  on  the  optical  properties  of  the  water. 

OBJECTIVES 

Optically,  phytoplankton  reveal  their  presence  through  their  influence  on  the  inherent  optical 
properties  (IOP’s)  of  the  water.  The  main  effect  of  phytoplankton  is  to  increase  the  absorption  of  light 
by  virtue  of  the  strong  absorption  by  their  photosynthetic  (chlorophyll  a)  and  accessory  pigments.  A 
secondary  effect  is  to  increase  the  backscattering  coefficient  of  the  medium  in  a  manner  that  depends 
on  the  concentration  of  pigments.  Although  techniques  for  measuring  the  absorption  coefficient 
directly  (e.g.,  in-situ  AC9  measurements  or  in-vitro  filterpad  absorption)  have  become  accepted  by  the 
scientific  community,  laboratory  techniques  for  measuring  backscattering  are  tedious  and  subject  to 
error,  and  in-situ  techniques  for  backscattering  are  in  their  infancy.  Thus,  there  has  been  considerable 
effort  devoted  toward  indirectly  inferring  these  IOP’s  by  virtue  of  their  influence  on  the  apparent 
optical  properties  (AOP’s),  e.g.,  the  diffuse  reflectance  of  the  water  (the  color  of  the  water)  or  the 
downwelling  irradiance  attenuation  coefficient.  These  AOP’s  are  perhaps  the  most  frequently 
measured  quantities  in  hydrologic  optics.  Clearly,  interpretation  of  such  observations  requires  a 
detailed  understanding  of  the  influence  of  phytoplankton  on  the  IOP’s,  and  their  link  to  the  AOP’s. 

The  IOP  <-»  AOP  link  forms  the  primary  focus  of  the  present  research.  In  particular,  our  research  is 
centered  on  deriving  the  IOP’s  from  measurements  of  the  AOP’s.  This  is  an  example  of  the  inverse 
problem  of  radiative  transfer.  It  is  important  in  that  IOP’s  determined  from  AOP’s  are,  by  definition, 
sampled  at  a  scale  appropriate  for  radiative  transfer,  and  for  remote  sensing.  Also,  the  retrieved  IOP’s 
possess  the  attribute  that  when  combined  with  the  radiative  transfer  equation,  they  reproduce  the 
measured  AOP’s. 

In  a  series  of  papers,  we  have  examined  the  inversion  of  apparent  optical  properties  (AOPs),  the  up 
(Eu)  and  downwelling  (Ed)  irradiances  or  upwelling  radiance  (Lu)  and  downwelling  irradiance,  to 
obtain  the  inherent  optical  properties  (IOPs)  -  vertical  profiles  of  the  absorption  (a)  and  backscattering 
(bb)  coefficients.  The  first  paper  reported  development  of  inversions  for  homogeneous  media  with 
elastic  scattering  [Gordon  and  Boynton,  1997],  the  second  extended  the  development  to  elastically- 
scattering  vertically- stratified  water  bodies  [Gordon  and  Boynton,  1998],  the  third  treated  inversion  in 
the  presence  of  the  interfering  effects  of  Raman  scattering  into  the  spectral  band  of  interest  [Boynton 
and  Gordon,  2000],  and  the  fourth,  provided  an  improvement  to  the  Gordon  and  Boynton  [1998] 
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algorithm  for  clear  water  (i.e.,  waters  in  which  scattering  by  the  water  itself  makes  a  significant 
contribution  to  bi,(z),  where  z  is  depth)  [Boynton  and  Gordon,  2002].  For  a  complete  program  of  AOP 
inversion  in  natural  waters,  the  only  process  we  have  overlooked  is  fluorescence,  i.e.,  inelastic 
scattering  into  a  given  band  of  wavelengths  from  all  shorter  wavelengths.  In  this  report  we  describe 
our  progress  on  the  inversion  of  AOPs  to  obtain  IOPs  in  fluorescing  media. 

APPROACH 


Consider  a  fluorescing  medium  described  by  IOPs  c(A),  /?(©,  A),and  b/(Xe—>X),  respectively,  the  beam 
attenuation  coefficient  at  the  wavelength  of  interest  A,  the  volume  scattering  function  for  elastic 
scattering  through  an  angle  0  at  A,  and  the  fluorescence  function  describing  the  fluorescence  from  an 
excitation  wavelength  Xe  to  A.  Assuming  the  fluorescent  emission  is  isotropic,  fluorescence  adds  the 
following  source  (intensity  density)  to  the  radiative  transfer  equation  (RTE): 

—  \bf  (z,Xe  -»  X)Eq(z,X6)  dXe, 

^ Xe<k 

where  Eq  is  the  scalar  irradiance  at  the  excitation  wavelength.  Henceforth,  we  will  use  standard 
terminology  and  symbols  of  ocean  optics  throughout  this  report,  so  not  all  symbols  will  be  explicitly 
defined.  For  a  fluorescing  medium,  we  add  to  a(z,X)  and  bb(z,X),  the  fluorescence  excitation  function 
bj(z,Xe^X),  as  the  IOPs  that  we  desire  to  retrieve  from  the  AOPs.  A  new  approach  for  inversion  in  the 
presence  of  inelastic  processes  was  outlined  by  Gordon,  [2002].  It  provides  a  partial  spectral 
decomposition  of  bj(Xe—>X). 


We  divide  the  excitation  spectrum  into  a  large  number  (N)  of  spectral  intervals,  the  ith  being  denoted  as 
A Xi  centered  at  A,.  Then  Gershun's  law  becomes: 

dE\  (z,X  ;)  N' 

=  -a{Xj  )Eq(z,  Xj )  +  X  bf(Xi^>Xj)EQ(z,Xi)AXi, 
i  =  1 


dz 


where  j  =  1  to  N,  N'<  N  is  the  number  of  intervals  use  to  fill  the  interval  between  Xe  and  A,  and  E\  =  E(i 
-  Eu.  We  have  assumed  that  the  entire  medium  is  homogeneous.  Writing  bj(z,Xi—>Xj)AXi  as  B(i,j )  and 
a(z,Xj)  as  we  can  write  Gershun’s  law  in  matrix  form: 


where 


d  E 


r  EQ(zr  1) 

EQ(zr  2)  •• 

-  EQ(zrj)' 

B  = 

B(2J) 

•  Ei  = 

Ey(ZTj) 

,  and  Eq  = 

Eq(z2X) 

Eq(z2,2) 

Eq^z2' j) 

[Ei(zrJ)j 

tw) 

EQ(zr2)  - 

•  Eo(zrj\ 

If  both  the  scalar  and  vector  irradiances  are  measured,  the  solution  is  found  by  inverting  the  scalar 
irradiance  matrix  Eq,  i.e., 


B  =  E 


-1 

0 


d  Ej 

dz 
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To  proceed  further,  and  estimate  the  backscattering  coefficient  we  need  to  solve  the  inverse  problem 
completely.  To  do  this,  we  note  that  the  quantities  a{Aj)  and  bj(Ai—>Aj),  A  i<  Aj,  are  already  known,  so 
only  the  volume  scattering  function  fi(z..0,A)  =  b(z,A)  P(z,0,A),  where  P(z,0,A )  is  the  elastic  scattering 
phase  function  of  the  medium  is  unknown.  If  we  assume  a  phase  function  for  the  medium,  then  only 
b(A)  is  unknown.  Thus,  one  need  only  solve  the  RTE  as  a  function  of  b(A)  and  choose  the  value  that 
provides  the  best  fit  to  Ed(z,A ),  Eu(z,A),  and  EAz.A).  The  backscattering  coefficient,  is  then  determined 
given  the  phase  function.  We  note  that  in  general  with  this  procedure  the  resulting  b/,(A)  is  almost 
completely  independent  of  the  assumed  P(0,A),  e.g.,  see  Gordon  and  Boynton  [1997]. 

It  is  important  to  recognize  that  the  solution  provided  above  does  not  require  any  knowledge  of  the 
IOPs  for  A  i<  Aj,  i.e.,  for  any  wavelength  smaller  than  Aj.  However,  it  is  required  that  all  wavelengths  A 
i<  Aj  that  contribute  to  the  fluorescence  at  Aj  are  included  in  the  analysis,  i.e.,  E\  and  Eq  (actually  only 
Eq)  are  measured  at  all  wavelengths  that  could  contribute  to  the  fluorescence  at  the  wavelength  of 
interest. 


If  Eo(ziJ)  is  not  measured,  it  must  be  replaced  by  Ji{z[ ,  j)~  '  E\  ( zi ,  j) ,  where  Ji{z[ ,  j)  is  unknown. 
Thus,  we  need  a  method  of  estimating  /7(z; ,  j) .  This  was  developed  based  on  the  assumption  that  the 
angular  distribution  of  the  fluorescence  is  uniform  in  the  upper  and  lower  hemispheres.  It  enables  one 
to  obtain  a  good  estimate  of  /j(z[  ,  j)  from  the  total  light  field  near  the  surface  (assumed  initially  to  be 
totally  elastic)  and  at  depth.  As  our  procedure  is  iterative,  we  will  use  the  superscript  “k”  to  indicate 

the  iteration  number.  Then  given  Jt^\z.[,  j)  we  can  form  the  kth  approximation  to  Eq^  and  with  the 
(k) 

new  Eq  produced  at  each  pass  yielding  a  new  estimate  of  B,  i.e., 


WORK  COMPLETED 


BW=[E^]  l^~. 
u  dz 


We  have  developed  a  computer  code  to  perform  the  above  inversion,  and  tested  it  with  synthetic  data. 

RESULTS 


We  examined  a  case  with  three  wavelengths:  A\  -  460  nm,  Ai  =544  nm,  and  A3  =  666  nm.  It  is 
assumed  that  there  is  no  inelastic  scattering  into  A\  from  shorter  wavelengths,  there  is  fluorescence 
from  A\  into  A2,  and  fluorescence  from  A2  into  A3;  however,  there  is  no  fluorescence  from  A\  into  A3. 
The  IOPs  associated  with  absorption  and  scattering  by  water  (aw  and  bw)  and  constituents  (ap  and  bp) 
that  were  used  are  provided  in  Table  1.  Simulated  light  fields  were  generated  at  each  wavelength  using 
a  Henyey-Greenstein  phase  function  (g  =  0.9)  for  elastic  scattering.  The  inelastic  scattering 
coefficients  used  are  provided  in  Table  2  (under  the  column  labeled  “True”),  and  the  AX s  were  taken 
to  be  1  nm.  The  goal  is  to  retrieve  the  IOPs  at  (and  into)  A3,  given  measurements  of  the  irradiances  at 
A\,  A2,  and  A3.  The  algorithm  described  above  was  operated  to  recover  values  of  atA-,),  bf(A\  —^3), 
bj{A2—>A3),  and  bXAi).  The  results  are  presented  in  Table  2  under  the  column  labeled  “Eq  -  Retrieved,” 
when  Eq  and  E\  are  measured  at  A3,  and  “£)  -  Retrieved,”  when  only  E\  is  measured  at  A3.  The  results 
of  the  inversion  are  excellent,  the  only  blemish  being  the  negative  (albeit  small)  value  retrieved  for 
bj(Al->A3). 
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Table  1:  Absorption  and  elastic  scattering  coefficients  (m1)  for  the  algorithm  example. 


IOP 

^3 

aw 

9.80  x  10'3 

5.02  x  10'2 

4.22  x  10  1 

bw 

4.12  x  10‘3 

2.06  x  10'3 

9.28  x  10'4 

Clp 

3.62  x  10'2 

1.60  x  10'2 

2.14  x  10'2 

_ K _ 

1.76  x  10  1 

1.49  x  10'1 

1.22  x  10  1 

Table  2:  True  and  retrieved  quantities  for  the  algorithm  example. 


Quantity 

True  (m'1) 

E()  -  Retrieved  (m'1) 

Ei  -  Retrieved  (m"1) 

a(X 3), 

0.4532 

0.4534 

0.4540 

bf) li  —>^2) 

381  x  10'b 

— 

— 

bf(A\ 

0 

-7.23  x  10'6 

-0.28  x  10'6 

bfA 2—>Af) 

144  x  10'b 

157  x  10'b 

144  x  10'b 

bb(A l3) 

2.80  x  10 3 

2.79  x  10'3 

2.82  x  10'3 

IMPACT/APPLICATIONS 

The  results  presented  above  suggest  that  the  algorithm  shows  promise  for  the  inversion  of  fluorescent 
light  fields  in  natural  waters.  We  are  presently  trying  to  understand  and  refine  its  performance,  and  to 
ascertain  its  usefulness  in  the  marine  environment.  We  believe  that  our  irradiance  inversion  algorithms 
will  be  of  significant  utility  for  processing  existing  and  future  experimental  irradiance  profile  data  to 
estimate  the  absorption  and  backscattering  coefficients,  and  their  relationship  to  constituent 
concentrations,  for  use  in  ocean  color  remote  sensing  algorithms.  In  particular,  these  methods  can  be 
used  with  older  data  sets  for  which  the  full  suite  of  IOPs  were  not  available. 

TRANSITIONS 

RELATED  PROJECTS 

We  are  now  collaborating  with  A.  Morel  to  apply  our  irradiance  inversion  algorithm  to  analysis  of  the 
OLIPAC  data  set.  We  are  also  collaborating  with  R.  Leathers  (NRL)  in  comparing  retrieval  algorithms 
operating  on  irradiance  data  from  the  Gulf  of  California,  and  with  the  PML  group  applying  our 
algorithms  to  data  acquired  in  coccolithophore  blooms. 
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le-5  le-4  le-3  le-2  le-1  le+O  le+1 


Figure  1:  Scalar  irradiance  at  the  three  wavelengths  for  the  example  provided 
in  the  text.  There  is  fluorescence  from  460  to  544  nm  and  from  544  to  666  nm, 
but  none  from  460  to  666  nm.  The  columns  of  the  Eo  matrix  contains  these 
profiles  (at  three  depths ).  The  nearness  of  the  460  and  544  nm  profiles 
suggests  the  Eq  matrix  may  be  close  to  singular. 


Fluorescent  Contribution  to  Irradiances 
0.0  2.0e-5  4.0e-5  6.0e-5  8.0e-5  1.0e-4  1.2e-4 


Figure  2.  Comparison  between  the  true  and  retrieved 
profiles  of  the  fluorescence  contribution  to  Efz)  and  Eu(z) 
for  the  example  in  the  text.  The  agreement  is  excellent. 
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